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We study the stationary and transient behaviors of the 
microemulsion phase subjected to a shear flow. The system 
is described by a diffusion-convective equation which general- 
izes the usual Cahn-Hilliard equation. Non-linear terms are 
treated in a self-consistent approximation. Shear, first and 
second normal stresses are calculated as momenta of the struc- 
ture factor. Shear thinning is observed in stationary condi- 
tions. After a newtonian regime at small values of the shear 
rate, the excess viscosity decreases when the shear rate be- 
comes of the order of the inverse of the relaxation time of 
the system without flow. In transient regimes, when the flow 
is applied starting from a quiescent state, we find that the 
shear stress reaches a maximum before decreasing to a con- 
stant value. 
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I. INTRODUCTION 

The rheological behavior of complex fluids such as 
polymer solutions, polymer melts, emulsions is of con- 
siderable interest both in technology and in basic re- 
search [1]. While the behavior of the stress response 
to applied flows is of fundamental importance in many 
applications, it also reflects the existence of mesoscopic 
structures in the fluid and is intimately related to its con- 
stitution. For example, when a shear flow is applied to a 
polymer solution, the stress first reaches a maximum and 
then relaxes to a constant value [2]. This phenomenon, at 
small strain, is related to the entanglement of the poly- 
mer network which is distorted by the flow with a result- 
ing increase of the stress. At larger values of the strain, 
however, the disentanglement of the system is favoured 
and the stress is observed to decrease. In general, non 
monotonic relaxational properties of the stress are typ- 
ical of complex fluids which are also characterized, in 
stationary conditions, by non-newtonian behavior. The 
effective viscosity depends on the applied shear flow and 
different behaviors can be observed [1]. 

In this paper we consider the rheological behavior of 
the microemulsion phase both in stationary and transient 
conditions. In ternary self-assembling systems the surfac- 
tant forms interfaces between oil- and water- like domains. 
These interfaces, in the microemulsion phase, constitute 
an intertwined bicontinuous structure disordered on large 
scales but with mesoscopic order on distances of the or- 
der of 500 Angstroms [3]. The observed structure factor 
is given by 



Hi) 



1 

a + gq 2 + cq 4 



(1) 



which, for g < 0, has a maximum at q = [4-6]. In 

real space, this corresponds to the two-point correlation 
function 
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where £ plays the role of the usual correlation length in 
disordered phases, d is related to the size of coherent 
regions of oil or water domains, and typical values of the 
ratio d/£ are in the interval 2 — 4. 

The behavior of the stress in the microemulsion phase 
was first considered in [7] and then in [8] where also 
the two-time correlation functions were studied. Here 
we complete the analysis of the steady state of [8] and 
consider also the transient behavior. Our approach is 
based on the use of a continuum free-energy functional 
and is similar to that of Onuki and Kawasaki [9] , applied 
also to evaluate the effects of a shear flow on copolymer 
melts [10,11], on the disorder-lamellar transition [13,14], 
and on the phase separation of binary mixtures [15]. 

We consider a Cahn-Hilliard equation generalized by 
the presence of a convective term. Hydrodynamical ef- 
fects are neglected; moreover, the surfactant is assumed 
to relaxe faster than the other components of the mix- 
ture so that its degrees of freedom arc not cxplicitely con- 
sidered. Non-linear terms, which become relevant close 
to transition lines, will be treated self-consistently. A 
renormalization procedure is introduced and the system 
is studied in terms of the physical variables £, d of the 
case without flow. 

Our main result for the stationary regime is the be- 
havior of the constitutive curve. Shear thinning, which 
is the decreasing of the effective viscosity when the shear 
rate is increased, is observed in two different ranges of 
the shear rate. It first occurs at a value of the shear rate 
of the order of the inverse of the relaxation time of mi- 
croemulsions without shear. The morphological changes 
occurring when the shear rate is increased can be deduced 
by the patterns exhibited by the structure factor. In the 
transient behavior after the application of the flow, at 
sufficiently high shear rates, we observe a maximum in 
the shear stress followed by a relaxation to a constant 
value, analogously to what observed in other systems. 
We have also studied the behavior of the stress tensor 
when, starting from a stationary state with shear, the 
flow is switched off. 
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The paper is divided as follows. In Section 2 we specify 
the model and solve formally the dynamical equation for 
the structure factor. Results for the stationary regime 
and for the transients are respectively described in Sec- 
tions 3 and 4. Section 5 contains some conclusions. 



II. THE MODEL 

Our study of rheological properties of microemulsions 
is based on a functional Landau-Ginzburg approach with 
one scalar order parameter <p(x) representing the concen- 
tration difference between oil and water. We consider the 
hamiltonian 



H[<f>] = Jd 3 x^ 
+ c(V 2 0) 



a 2 (, 



(3) 



which has been largely demonstrated to well describe 
equilibrium properties of ternary mixtures [3] . Here we 
briefly discuss the properties of this hamiltonian relevant 
for the microemulsion phase. The expression (3) differs 
in the gradient terms from the usual Landau-Ginzburg 
hamiltonian used to study binary mixtures. A negative 
value of the function g(<j>) = go + .92<^ 2 (52 > 0) favors 
the appearing of interfaces. In particular the value of g 
can be related to the amount of surfactant present in the 
system. The term proportional to c > assures stability 
at large momenta and weights the curvature of interfaces. 

The presence of the quartic terms 4 and 4> 2 , 

which also in a disordered phase could have a role in the 
proximity of a transition line, makes impossible an ex- 
act determination of the two point correlation functions. 
However, following [7], it is possible to use a renormaliza- 
tion procedure based on a self-consistent approximation 
to find an expression for the equilibrium scattering func- 
tion and for the two-point correlation function in real 
space. For the first, defined <j>(k) the Fourier transform 
of <f>(x), it is found that 



S{k) =< <f>(k)<f>(-k) >- 
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where k = \k\, the renormalized parameters are given by 

a r = a 2 + XSq + g 2 S 2 (5) 
9t = 5o + 52^0 (6) 

and the loop integrals are defined as 
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(p = 0, 2) with A being a high momentum phenomcno- 
logical cut-off. The region g r < and 4a r c — g 2 > 0, 



as discussed in the introduction, can be identified with 
the microemulsion phase with the function S(k) having 

a peak at k = &m = \J %r- Moreover, the characteris- 
tic lenghts £ and d appearing in the real space two-point 
correlation function (2) are given by 
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Therefore equilibrium properties can be expressed in 
terms of renormalized parameters and, by equations (8) 
and (9), through the physical lenghts £ and d (once the 
parameters A, c, g 2 and the cut-off A are given). 

The dynamics of the order parameter in presence of 
convective motion [16] is described by the equation: 



af + «-v* = rv Hj-, 
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where H is the hamiltonian of (3) . The velocity field v is 
a planar Couette shear flow: 
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where 7 is the shear rate and e x the unit vector in the 
flow direction; rj^ is a white gaussian noise representing 
thermal fluctuations with momenta given by 

< V<t>(x,t) > = (12) 
< 770 (x, t)»fo(f , t') > = -2TTV 2 S 3 (x - x')5{t - t') (13) 

(< • • • > means the ensemble average) as required by the 
Fluctuation-Dissipation theorem which holds in absence 
of flow. The functional derivative 5H/5(f> represents the 
difference in chemical potentials between oil and water; 
r is a mobility coefficient and T is the temperature of 
the heat bath. By assuming equation (10) as the evo- 
lution equation for 4>, we are neglecting hydrodynamics 
fluctuations as well as the motion of the surfactant. 

We will study the evolution equation for the dynamical 
structure factor 



S(k,t) =< (j>(k,t)<p(-k,t) > 



(14) 



in the same self-consistent approximation used in equilib- 
rium to write Eqs. (5,6). The convection-diffusion equa- 
tion can be formally linearized as [17] 

+ v- V(f> - r V 2 {(a 2 + A5b(i) + 92S 2 (t)) c\> 

- (g o + g 2 S o (t))V 2 ^ + cA 2 0}+^ (15) 

where the quantities S p (t) are given by expressions ana- 
logue to those of Eq. (7) but now with S(k, t) of Eq. (14) 
self-consistently calculated with Eq. (15). A standard 
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procedure gives from Eq. (15) the dynamical equation 
for S(k,t): 



In addition, the excess viscosity is defined as: 



d 
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where K R {k) = a r +g r k 2 + ck A is the renormalized vertex 
function and the parameters a r and g r can be obtained 
as in Eqs. (5,6) using S p (t). 

A formal solution of Eq. (16) may be obtained by the 
method of characteristics: 

S(k,t) = A Q (lC(t))l 1 (t) + 2TTl 2 (t) (17) 

where we have defined the functions 

2^ = e -2r f* dsK 2 (s)[a r + ar K 2 (s)+cK l (s)] ^ 

2 2 (t) = f du/C 2 (w)Ji(w) (19) 
Jo 

and IC(u) = k + jk x ue y ; Ao(fc) is the structure factor at 
the initial time t = 0. Since the quantities a r and g r con- 
tains the momenta of S(k, t) 1 Equation (17) is actually a 
nonlinear integral equation for S(k, t). This equation can 
be solved numerically at each time by iterative methods. 

Our results will first concern steady state properties. 
The stationary solution can be readily obtained from the 
t — ► +oo limit of Eq. (17), observing that in this limit 
the first term of the solution tends to zero (except for the 
k = mode). Therefore we write the stationary structure 
factor as: 



S(fc;7)oo = 2Trj 2 (oo) 
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where 

J 2 (oo) = I™ dzK,\z)e- 2T /;« 2 («)K+^ 2 W+^( S )] 
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We will also study transient behaviors with the fluid 
evolving from a quiescent state towards the stationary 
state with shear, or with the system relaxing, after in- 
terruption of the flow, from the sheared stationary state 
into the quiescent state. For the latter case we use the 
solution of equation (16) with 7 = 0: 

S(k',^f',t) re i ax = S(k; r ^) 00 e R ^ ' 

+ _l_( 1 _ e -^^T) (21) 

Finally, once the structure factor is known, we may 
evaluate the stresses which can be obtained as momenta 
of the structure factor [7]. The shear, first and second 
normal stresses are respectively given by 
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which represents the contribution of interfaces to the full 
viscosity of the fluid (that is, evaluating the viscosity of 
the fluid by means of equation (25) we are neglecting the 
hydrodynamical contribution to the viscosity itself). 



III. STATIONARY REGIME 

In this section we present results for the steady states 
reached under the action of the shear flow with the struc- 
ture factor given by Eq. (20). We have studied numeri- 
cally this expression for several values of £, d and 7. The 
other parameters have been fixed as g 2 = l,c = I, \ = 
0.5, A = 3. 

The effects of the flow on the structure factor can be 
seen in Fig. 1 where the projections on the planes k y = 
and k z — are shown for different 7 and £ = 2, d = 6. 
Similar results have been obtained for other choices of 
£ and d [8]. (At k x — the shape of the structure fac- 
tor is the same of the case without flow, see Eq. (16)). 
At 7 = 0.5 the structure factor remains almost isotropic 
and its pattern for each cartesian plane is close to that 
of a circular volcano. The patterns are progressively dis- 
torted when the shear rate is increased. On the plane 
k z = 0, at 7 = 2, the edge of the volcano has assumed an 
elliptical shape and four peaks are visible. These peaks 
initially appear on the coordinate axes; then, when 7 
is increased, the ones located at k x ~ become com- 
paratively more important while the two others rotate 
clockwise and decrease their amplitude linearly with 7 
until they disappear. Indeed, in the limit 7 — > 00, since 
terms proportional to powers of jk x damp the exponen- 
tial term on the r.h.s. of Eq. (20), only the maxima of 
C{k) with k x = and k y = ikM survive. On the other 
plane k y = two peaks at k x — 0, k z = ±kM, are also 
observed to become sharper and sharper as 7 is increased. 

The above results can be related to the orientation 
of the interfaces in the mixture, as also observed in [8]. 
A peak of C(k) defines a characteristic length propor- 
tional to the inverse of its position and, since the system 
is not isotropic, to each maximum one associates three 
lengths, one for each space direction. Due to the sym- 
metry k — > —k only the peaks not related by reflection 
around the origin can be considered. At very large shear 
rate the existence of a single couple of maxima at k x — 
signals that interfaces are preferentially aligned along the 
flow with symmetry recovered in the transverse directions 
and the characteristic lengths being the same as without 
shear. For .intermediate values of 7 the additional peaks 
at {k x ,k y ,k z ) reveal the presence of interfaces oriented 
with an angle a = a,rcta,n(—k x /k y ) with respect to the 
flow, besides those aligned along the x direction. These 
additional peaks are better seen in a region of parameters 
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closer to the microemulsion-laniellar transition line cor- 
responding to a larger value of £ [8] . As 7 is increased the 
tilt angle a and the relative abundance of lamellae orien- 
teted at this angle diminish as suggested by the behavior 
of the maxima with k x ^ previously discussed. 

The behavior of the stress tensor as a function of the 
shear rate is reported in Figs. 2-4. At small values of 
7 the shear stress is a linear function of 7 so that the 
viscosity is constant and the fluid is ncwtonian. Shear 
thinning occurs for 7 between 1 and 2, when the slope 
of the curve of the stress changes significantly. At this 
point the original volcano shape of the structure factor 
has also appreciably changed. In terms of interfaces we 
expect that, when shear thinning is observed, the bicon- 
tinuos network of interfaces which is only distorted in the 
newtonian regime, is affected by many ruptures with a 
significant dicrease of connectivity. We can say that a 
strong shear regime is entered. Indeed, we can compare 
the shear temporal scale r y~ 1 with the relaxation time t m 
of microemulsions in equilibrium. The relaxation time of 
a mode with wavevector k is given by 



r(k) = 



1 



Tk 2 (a r +g r k 2 +ck 4 ) 



(26) 



Following [7] we choose k — ku corresponding to the 
peak of microemulsions, so that 
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The corresponding Deborah number is given by 



TS 



(27) 
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where Tg = I/7. If we take 7 = 2 we get De ~ 1.9 
for the case of Fig. 2. This indicates that shear thinning 
becomes evident when the shear rate is of the order of 
the inverse of typical structural times of the system with- 
out shear. We checked for other values of £, d that the 
Deborah number at shear thinning is always of order 1. 

Shear thinning is observed again at 7 ~ 10 3 when, as 
we have seen, peaks with k x = largely prevail. We 
expect that at these values of the shear rate the original 
bicontinuous interface network has changed significantly 
its topology becoming more similar to a stack of lamellae. 
At very large 7 the excess viscosity is found to decrease 
as l/7~ s with s = 1.87 which is close to the analytical 
limit s = 2 [7]. In this final stress regime the lamellae 
are expected to become more and more aligned with the 
flow with fluctuations very inhibited. We observe that 
the shear stress corresponding to a completely ordered 
lamellar phase is zero. 

The other stress components N\ , N2 behave similarly. 
At small 7, Ni,N 2 ~ 7 2 while they decrease as 7 _1 when 
7 — > 00 (see Figs. 3-4). 



IV. TRANSIENTS 

We have studied the evolution of the system under 
the action of the shear flow from the initial equilibrium 
configuration of Eq.(4) towards the steady state of the 
previous section, as described by Eq. (17). The behavior 
of the stress components for different values of 7 is shown 
in Figs. 5-7. When 7 is large enough that De > 1, a 
non monotonic behavior of the stress is observed with 
&xy, Ni, IJV2I exhibiting a maximum before relaxing to a 
constant value. A similar behavior has been measured in 
polymer solutions [2]. In our case we can think that at 
initial times the surfactant interfaces are stretched by the 
flow with a consequent increase of the stress. When the 
maximum of the stress is reached, the interface structure 
starts to be broken and the stress relaxes to a lower value. 
The temporal evolution of the structure factor for 7 = 
100 is shown in Fig. 8. The largest observed distortion 
corresponds to the maximum of the stress. In the case 
7 = 2, when the relaxation of the stress is monotonic, a 
prolate pattern like that in the middle of Fig. 8 at k z = 
is not observed. 

We have also considered the opposite situation with the 
system, initially in a stationary state with shear, evolving 
without flow as described by Eq. (21). In this case the 
behavior is exponentially monotonic after an initial faster 
decay, as it can be seen in Fig. 9. The time constant r 
of the exponential part of the relaxation decreases with 
7, as shown in Fig. 10. 



V. CONCLUSIONS 

We have used a generalized Cahn-Hilliard equation 
with a convective term to study the rhcological behav- 
ior of the microemulsion phase. The steady state con- 
stitutive curve shows shear thinning first occurring at a 
shear rate of the order of the inverse of the equilibrium 
relaxation time. We have also obtained for the first time 
analytical expressions for the temporal behavior of the 
structure factor. From this we derive a non monotonic 
evolution of the stress. This is similar to what is ob- 
served in other systems that relaxe into the steady state 
with a shear flow. We hope that these predictions are 
useful for future experiments. From the theoretical point 
of view this analysis can be completed studying how the 
equilibrium phase diagram, including the disorder and 
the Lifshitz lines [3-6], is changed by the presence of 
the flow. Moreover, hydrodynamic fluctuations should 
be taken into account for a full description of the sys- 
tem. 
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FIGURES 



FIG. 1. Projections of the structure factor in the stationary state on the planes k y = (right column) and k z — (left 
column) for £ = 2, d — 6. The shear rate is, from the top to the bottom, 7 = 0.5, 7 = 2, and 7 = 100; k x ,k y , k z vary between 
—3 and 3 in adimensional units. 
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FIG. 2. 



The stationary shear stress a xy versus the shear rate 7 for £ = 2 and d 




FIG. 3. The first normal stress Ni for £ = 2 and d — 6. 
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FIG. 4. Absolute value of the second normal stress N2 for £ = 2 and d — 6. 
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FIG. 5. Time evolution of the stress a xy for various values of the shear rate 7 with £ = 2 and d = 6. Results are shown for 
7 = 0.5 (top- left), 7 = 2 (top-right), 7 = 20 (bottom- left) and 7 = 100 (bottom-right). 
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FIG. 6. Time evolution of the first normal stress Ni; the parameters are the same of Fig. 5. 
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FIG. 7. Time evolution of the second normal stress N2 with parameters as in Fig. 5. 
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FIG. 8. Time evolution of the structure factor for £ = 2, d — 6, and 7 = 100. The projections on the planes k z = (left 
column) and k y — (right column) are shown respectively for t = (top), f = 2.5 x 10~ 2 (middle), when the maximum of a xy 
is reached, and t = 2.0x 10 _1 (bottom). 
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FIG. 9. Relaxation of the stress a xy for £ = 2, d — 6, 7 = 2 (top-left), 7 = 20 (top-right) and 7 = 100 (bottom). 
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FIG. 10. Relaxation time r as a function of the shear rate 7. 
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